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Abstract 

Models of folding of a triangular lattice embedded in a discrete space are studied as simple 
models of the crumpling transition of fixed-connectivity membranes. Both the case of planar folding 
and three-dimensional folding on a face-centered-cubic lattice are treated. The 3d-folding problem 
corresponds to a 96- vertex model and exhibits a first-order folding transition from a crumpled phase 
to a completely flat phase as the bending rigidity increases. 

1 Introduction 

The statistical mechanics of polymers, essentially one-dimensional objects, has proven to be a rich and 
fascinating field. |^] The success of physical methods applied to polymers arises from universality — 
many of the large-length-scale properties of polymers are independent of microscopic details such as the 
chemical identity of the monomers. The generalization of polymer statistical mechanics to membranes, 
two-dimensional surfaces fluctuating in some embedding space, has proven to be even richer and is still 
under active development. In contrast to polymers there are different universality classes of membranes 
distinguished by their long-range order. These classes are the analogues of the well-known crystalline, 
hexatic and fluid phases of strictly two-dimensional systems. ^] This diversity arises from the richer 
geometry of surfaces as compared to lines and the resultant enlarged space of possible symmetries. 

The closest membrane analogue to a polymer is a 2D fishnet-like mesh of monomers (or vertices) 
with fixed connectivity. The bonds (links) are assumed to be unbreakable. The vertices themselves live 
in R d , with a physical membrane corresponding to the case d = 3. Such a membrane is called a 
crystalline or polymerized membrane by virtue of its intrinsic crystalline order. In this talk we will 
deal only with phantom (non-self-avoiding) membranes. In general the Hamiltonian for a polymerized 
membrane will have both in-plane elastic (strain) contributions and out-of-plane bending contributions, 
since it can both support shear and fluctuate in the embedding space. In Monge gauge the most general 
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effective continuum Hamiltonian at long wavelength takes the form 

H{h,u) = % J d 2 a{V 2 hf + \ J d 2 a(2»u 2 af3 + \u 2 J, (1) 

where h is the height above a reference plane with intrinsic coordinates ((71,(72), u are the phonon 
modes, u a p is the strain tensor, k is the bending rigidity and \i and A are the bare Lame coefficients. 
The strain tensor u a p measures the difference between the induced metric and a flat reference metric 
determined by the equilibrium configuration of the membrane at rest. 

The distinctive physics of polymerized membranes arises from an unavoidable nonlinear cou- 
pling between the in-plane strain modes and the out-of-plane height fluctuations. Integrating out the 
quadratic phonon fluctuations by linearizing the strain tensor one finds an effective long-range interac- 
tion between Gaussian curvature excitations which tends to stiffen the surface on long length scales. || [| 
In other words, the effect of undulations on small length scales is to increase the bending rigidity on 
longer length scales. This effect is easily demonstrated — an ordinary piece of paper becomes much 
stiffcr with respect to bending deformations after it is crumpled and then opened up. More technically, 
one finds that the bending rigidity at sufficiently long wavelengths is momentum dependent n(q) ~ q v , 
with rj positive. This scaling behaviour leads to a stable flat phase for polymerized membranes,^, [| 
with remarkable properties controlled by a infrared stable fixed point. || For small bending rigid- 
ity (high-temperature), on the other hand, polymerized membranes are crumpled. Q At some critical 
bending rigidity (or equivalently critical temperature) there should, therefore, be a crumpling transition 
from a flat to a crumpled phase. 

The flat phase of polymerized membranes is characterized by long-range orientational order 
in the surface normals. Since long-range order is unusual in 2D systems it is worthwhile to explore 
as many avenues as possible to understand its exact nature and origin. Aronovitz and Lubensky || 
analyzed the flat phase of D-dimensional elastic solids embedded in R d , and subject to bending energy, 
using an e expansion about the upper critical (manifold) dimension D uc = 4, with fixed codimension 
d c = d — D. They determined the renormalization group flows in the space of dimensionless couplings 
fi = /j,l e /K 2 and A = Xl e /n 2 , where I is the renormalization length scale. They discovered a globally 
attractive infrared fixed point describing a flat phase and determined the scaling dimension rj and the 
corresponding exponent r\ u (fi ~ A ~ q nu ) for the infinite renormalization of the Lame coefficients. 

A revealing extreme limit of polymerized membranes was studied by David and Guitter. [fl0| 
This is the limit of infinite Lame coefficients in the Hamiltonian (0). Since the elastic term scales like 
q 2 , in momentum space, as compared to q 4 for the bending term, this limit may be regarded as a means 
of determining the dominant infrared behaviour of (Q) . In this "stretchless" limit the strain tensor u a p 
must vanish and the model is constrained, very much in analogy to a non-linear sigma model. The 
/3-function for the suitably rescaled inverse bending rigidity a — d/n may then be computed in the 
large-<i limit, yielding 

„, . da 2 , 1 const. , . . 

W=«^=d a -fe + — (2) 
For d = 00 there is no stable fixed point and the membrane is always crumpled. Eq. ^ reveals, 
to next order in 1/d, an ultraviolet stable fixed point at a — Hw/d, corresponding to a continuous 
crumpling transition. The critical exponents associated with the "flat" fixed point, discussed above, 
may also be determined in the large-<i expansion. (7]] 



2 Folding 

A natural lattice formulation of the infinite-elastic-constant limit is the statistical mechanics of folding of 
a regular triangular lattice. A folding in R d of the regular triangular lattice is a mapping which assigns 
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Figure 1: The oriented triangular lattice: triangles pointing up (resp. down) are oriented counterclock- 
wise (resp. clockwise). The three tangent vectors tj, i = 1,2, 3, have a vanishing sum in the embedding 
space. 

to each vertex v of the triangular lattice a position X„ in the <i-dimensional embedding space R d , with 
the "metric" constraint that the Euclidean distance |X„ 2 — X„j | in R d between nearest neighbours v\ 
and V2 on the lattice is always unity. Under such a mapping, each elementary triangle of the lattice is 
mapped onto an equilateral triangle in R d . In general, two adjacent triangles form some angle in R d , 
i.e. links serve as hinges between triangles and may be (partially) folded. Folding is best described 
in terms of tangent vectors, which are link variables defined as follows: we first orient the links of the 
lattice as in Fig|l], with triangles pointing up (resp. down) oriented counterclockwise (resp. clockwise), 
and define the tangent vector between two neighbours v\ and v 2 as the vector t VltV2 — X„ 2 — ~K Vl if 
the arrow points from v\ to v 2 . The metric constraint states that all tangent vectors have unit length. 
With our choice of orientation, moreover, the three tangent vectors t^, i = 1,2,3, around each face of 
the lattice must have vanishing sum. This is the basic folding rule: 

ti + t 2 + t 3 = 0. (3) 

Up to a global translation in R d , a folding is therefore a configuration of unit tangent vectors defined 
on the links of the lattice, obeying the folding rule (||) around each triangle. 

Viewed in terms of normals membrane models resemble non-linear sigma models. In the lattice 
version the correspondence is with Heisenberg spin models. There is one key difference — apparent in 
both the continuum and lattice models. Normal vectors are not arbitrary unit vectors in the embedding 
space — they are constrained to be normal to the underlying 2D manifold. This means they are actually 
Grassmannian sigma models in the continuum |l2] | and constrained Heisenberg models on the lattice. 
These constraints play a crucial role in stabilizing an ordered phase. 

3 Planar Folding 

The two-dimensional or planar (d — 2) folding problem of the triangular lattice was first introduced 
and studied numerically by Kantor and Jaric.(]l| It is clear that the only remnant of the normal vector 
degree of freedom is a Z 2 spin corresponding to orientation (up or down). As an Ising spin system 
the planar folding model is constrained (as noted above) and is described by an 11-vertex model. In 
folding terms these 11 vertices describe the distinct folded states of a single hexagon of the original 
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Figure 2: The phase diagram for planar folding in the (K,h) plane, where K is the discrete bending 
rigidity and h is an external magnetic field. Three first order lines h = h c (K),— h c (K) (K < K c ) and 
h = (K > K c ) separate the three phases M = 0,±1 and meet at the triple point (K c ,0). 



triangular lattice. In terms of the total magnetization M the analogue of the crumpling transition for 
this model would be a spin-ordering transition from an unmagnetized phase (M = 0) to a magnetized 
phase (M ^ 0). Analytic progress on planar folding was made subsequently. jHj For planar folding 
it is easy to check that, up to a global rotation in the embedding plane, all the link variables are 
forced to take their values among a fixed set of three unit vectors with vanishing sum. This permits a 
reformulation of the pure planar folding problem (with no bending rigidity) as that of the 3-colouring 
of the links of the triangular lattice: calling the three fixed vectors blue, white and red, the folding 
rule (|J) translates into the constraint that the three colours around each triangle be distinct. This 
3-colouring problem was solved by Baxter ]l5|] with Bethe Ansatz and transfer matrix techniques. His 
result for the thermodynamic partition function measures the number of distinct folded configurations 
Z2d oc q^t f° r a lattice with Na triangles, in the limit of large Na- This gives the folding entropy per 
triangle s 2 d = log(g 2 d), with Q 

Q2 d = ^r(l/3) 3 / 2 = 1.20872... (4) 

The 2d folding problem has also been studied in the presence of bending rigidity, which associates an 
energy to each folded link, and with a magnetic field coupled to the sum of the normal vectors to the 
lattice. Ji6| The model was found to undergo a first order folding transition. At zero rigidity and zero 
magnetic field, the lattice is in an entropic folded phase (M = 0). At large enough rigidity and/or 
magnetic field, the lattice becomes totally unfolded (M — ±1). The phase diagram is shown in Fig. 



4 Three-Dimensional Folding 

The existence of a unfolded (magnetized) phase and a first-order folding transition is fascinating, but 
one would like to know how sensitive these features are to the discretization of the space of local normals. 
After all, the lower critical dimension for systems with discrete symmetry is typically less than that 
for systems with continuous symmetry. It is conceivable that the transition disappears altogether 
if we enlarge the embedding space or that the order of the transition changes. In particular there 
is no rigorous analytic prediction for the order of the crumpling transition for physical polymerized 
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Figure 3: The specific heat versus bending rigidity from the Monte Carlo simulation of Ref. 24. The 
critical exponent a = 0.4(1). 



membranes. Paczuski, Kardar and Nelson [[L7j found that critical fluctuations about the mean field 
solution, in an e = 4 — D expansion, drive the transition first order for embedding space dimension 
d < 219. While the early numerical simulations could not rule out a weak first-order transition, |l^] later 
more extensive simulations clearly indicate a continuous transition. j^, |2(| [2l], |||, |24| The specific 
heat plot from Ref. 24 is shown in Fig. To gain some insight into the effect of a finer discretization of 
the embedding space we have generalized planar folding to folding on a 3d Face Centred Cubic (FCC) 
lattice. [H) 

In the general 3-dimensional folding problem, the local folding constraint (|J) imposes only that 
the three tangent vectors around each face be in the same plane and have relative angles of 2n/3. This, 
however, does not impose any constraint on the relative positions of the two planes corresponding to 
two adjacent faces, which may form some arbitrary continuous angle. As opposed to the 2d case, this 
then leads to a problem with continuous degrees of freedom. 

To define a discrete model of folding in 3c?, one must allow only a finite number of relative angles 
between adjacent faces. More specifically, one may also impose that the link variables themselves take 
their values among a finite set of tangent vectors, now in R 3 . For symmetry reasons, we took this 
set of tangent vectors to be the (oriented) edges of a regular solid of R 3 , made of equilateral triangles 
only. There are only three regular solids in R 3 made of equilateral triangles: the tetrahedron, the 
octahedron and the icosahedron. The edges of the tetrahedron (resp. icosahedron), however, cannot be 
consistently oriented such that the corresponding tangent vectors satisfy (|j) around each face. This is 
because each vertex is surrounded by an odd number 3 (resp. 5) of triangles. There is no such problem 
for the octahedron, as shown in Fig. [|. The 12 links of the octahedron are oriented consistently to form 
8 triplets of tangent vectors with vanishing sum, corresponding to the 8 faces of the octahedron. One 
may therefore consider the restricted 3c? "octahedral folding" problem, where the tangent vectors are 
chosen from the set of 12 edge vectors of a regular oriented octahedron. In the folding process, the 
folding rule (|J) imposes that the three links of a given face on the original triangular lattice are mapped 
onto one of the 8 triplets of tangent vectors above. For a given triplet, the triangle can still be in 3! 
states corresponding to the 3! permutations of the three edges. Each triangle can therefore be in one 
of 48 = 8 x 6 states. 
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Figure 4: The oriented octahedron: the edges around each face form triplets of tangent vectors with 
vanishing sum. The four normal vectors n^, i = 1,2, 3, 4, are represented on the corresponding outward 
oriented faces. 

The 8 faces of the octahedron can be labelled as follows: we consider for each face its normal 
vector, pointing outward or inward according to the orientation of its tangent vectors on the octahedron 
(see Fig. |]). On the octahedron there are four alternating outward and inward oriented faces. The 
normal vectors to opposite faces are equal. This defines a set of four vectors ni, 112, 113 and 114, which 
furthermore satisfy the sum rule ni + ri2 + 113 + 114 = 0. Each face is labelled by its orientation (outward 
or inward) and its normal vector (1, 2, 3 or 4). 

The 12 oriented edge vectors of the octahedron are actually identical to the 6 edge vectors of 
a tetrahedron, now taken with both possible orientations. The four normal vectors above are also the 
normals to this tetrahedron. For each folding map, the image of the folded lattice in R 3 lies therefore 
on a 3d FCC lattice, which consists of a packing of space by octahedra complemented by tetrahedra, 
as shown in Fig. |5|. In this respect, the "octahedral folding" problem simply corresponds to discretizing 
the embedding space as an FCC lattice. 

4.1 96-vertex Model 

When stated in terms of tangent vectors, the 3d "octahedral folding" problem involves three types of 
constraints: face, link and vertex constraints. The first constraint, around each face, imposes that the 
three tangent vectors of a given triangle form one of the 8x6 (ordered) triplets with vanishing sum. 
The second constraint, on each link, arises because two adjacent triangles share a common tangent 
vector. Given the state of one triangle, any adjacent triangle has one of its tangent vectors already 
fixed and thus is left with only 4 = 48/12 possible states. They correspond simply to the four values 
for the relative angle between two neighbouring triangles, i.e. the angle between the normal vectors, 
depicted in Fig. |[ These four values are (no fold: the triangles are side by side), 180° (complete fold: 
the triangles are on top of each other), cos~ 1 (— 1/3) ~ 109°28' (fold with acute angle: the two triangles 
lie on the same tetrahedron) and cos _1 (l/3) ~ 70°32' (fold with obtuse angle: the triangles lie on the 
same octahedron). Finally, there is a third constraint on the six successive folds around each vertex 
of the lattice: after making one loop, the same tangent vector must be recovered. Since the "metric 
constraint" is local, there are actually no constraints other than these three (face, link and vertex) 
constraints. 

In the study of the 2d folding (3-colouring) problem, the face and link constraints are taken 
into account by going to Z2 spin variables defined on the faces of the lattice. Ordering the colours 
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Figure 5: The FCC lattice viewed as a packing of 3d space with octahedra and tetrahedra. 



acute 




Figure 6: The four possible folding angles between two adjacent triangles. The neighbour of the dark 
triangle may (i) be itself on top of the dark triangle (complete fold), (ii) occupy the symmetric position 
in the same plane (no fold), (iii) lie on the same octahedron (i.e. form an obtuse angle) or (iv) lie on 
the same tetrahedron (i.e. form an acute angle). 
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Figure 7: The transition from a link to a subsequent link (k,l) is described by the two Z2 face 
variables z and a. 

cyclically, the spin is +1 (resp. —1) if the colour increases (resp. decreases) from one link to the neigh- 
bouring one on the triangle, oriented counterclockwise. In this language, the actual folds take place 
exactly on the domain walls of the spin variable. Instead of having a Z3 colour variable per link, one 
is left with a Z2 spin variable per triangle. The vertex constraint translates into a constraint on the 
six spins &\, ...,06 around each vertex of the lattice, namely that °i = mod 3. This leads to 

22 possible local spin configurations around each vertex, or equivalently, after removing the global Z2 
degeneracy of reversal of all spins, to an 11-vertex model on the lattice. [|l4| 

One may proceed in the same way for the 3d "octahedral folding" and account for the face and 
link constraints by expressing folded configurations in terms of two Z2 variables on the triangles. These 
variables will indicate the relative states of successive links around the face. One may then count the 
number of allowed hexagonal configurations around a vertex: this leads to a 96-vertex model. These 
vertices and the corresponding rules on the Z2 variables will be given below. 

Let us label the 12 edges of the octahedron as follows: each edge is shared by two adjacent faces, 
one outward and one inward oriented (see Fig. We label the edges by the indices 1 < i 7^ j < 4, 

when the normal vector for the outward face is and the one for the inward face is rij. There are 12 
such couples Consider now an elementary triangle of the lattice. Starting from one of its links 

the subsequent link (k,l) counterclockwise must share a face with on the octahedron. This 
leads to the 4 following possibilities, labelled by the Z2 face variables z, a £ {±1}: 

z = +1 : -> 1 + j, €iji = -a = ±l 

z= -1 : -> (k,j), k^i, €i jk = +a = ±l, 

where djk — X)z tijki is defined in terms of the totally antisymmetric tensor eijki, equal to the signature 
of the permutation (ijkl) of (1234). The value z = +1 (resp. z = —1) indicates that the two tangent 
vectors share an outward oriented (resp. inward oriented) face on the octahedron. The spin variable 
a takes the value +1 (resp. —1) if (k,l) follows (resp. precedes) on their common (oriented) face 
of the octahedron. The variable a also indicates whether the normal vector to the triangle (in the 
embedding space R 3 ) is parallel (a = +1) or antiparallel (a — — 1) to the corresponding normal vector 
of the octahedron. 

Considering now two neighbouring triangles, the 4 possible relative values z-ijz\ and Gija\ 
indicate which type of fold they form. The domain walls for the z variable are the location of the folds 
which are either acute or obtuse, whereas those for the a variable are the location of the folds which 
are either complete or obtuse. The superposition of these two types of domain walls fixes the folding 
state of all the links, specifying the folding state of the lattice up to a global orientation. 

The use of z and a variables instead of the 12 variables incorporates the face and link 

constraints. As in the 2d case, the vertex constraint is more subtle. Nevertheless, one can count the 
number of possible configurations around a vertex satisfying this constraint, i.e. the number of possible 
folded states of an elementary hexagon. Indeed the mapping (||) may be represented by a 12 x 12 
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Figure 8: The six Zi and cr, variables around a given vertex, and the colours of the interior links. 

connectivity matrix M{i,j),(k,l) with i ^ j and k =/= I: 

M(i, 3 ),(k,i) = $ik + Sji - 2SikSji . (6) 

This matrix acts as a transfer matrix between two successive internal links of the hexagon. The number 
of configurations of a hexagon is simply given by: Tr(M 6 ) = 4608, where the trace guarantees that the 
same link variable is recovered after one loop. These 4608 configurations count as distinct all the foldings 
which are related by a global change of orientation of the hexagon in embedding space. The order of the 
resultant degeneracy is 48, corresponding to 12 choices for the first tangent on the octahedron times 4 
for the choice of the second from among its 4 neighbours (this latter choice corresponds to the 4 choices 
of the z and a variables on the corresponding triangle). This leaves us with 4608/48 = 96 distinct 
configurations. 

4.2 Vertex Rules 

Consider an elementary hexagon of the triangular lattice. As shown in Fig. || each of its 6 triangles is 
labelled by z and a variables. Let us also assign a colour c, to each link: Cj = Cq + Y]j—i o~j mod 3. 
Geometrically each colour corresponds to one of the three orthogonal planes of the target octahedron 
(comprising four links). The vertex constraint arises from the requirement that six applications of the 
link mapping of Eq. (|^) be the identity. By relating the link mappings to the elements of the tetrahedral 
group A\ one finds two folding rules. [^5| The first folding rule: 

6 

°i = mod 3 (7) 

for the six spins around the central vertex of the hexagon is identical to the planar folding rule. In 
contrast with the 2c? situation, the restriction (0) is not the only constraint. There is a second folding 
rule 

6 1 

^(l — z i z i+i)d( c ii c m °d 3) = mod 2 : c = 1, 2 . (8) 

8=1 

The two folding rules (Q) and ^ characterise the vertex constraint entirely. One finds 384 = 96 x 4 
vertex configurations (there is a 4-fold global degeneracy under reversal of z or a). The 96 folding 
vertices are displayed in Fig. |9| with the following conventions: no line corresponds to no fold; a thick 
line corresponds to a complete folding (180°, flip of a only); a thin line corresponds to a fold with 
obtuse angle (cos _1 (l/3) ~ 70°32' between normal vectors, flip of both a and z) and finally a dashed 
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Figure 9: The 96 vertices satisfying the two folding rules [7]and[|: no line corresponds to no fold, a thick 
line corresponds to a complete fold, a thin line corresponds to a fold with obtuse angle and a dashed 
line corresponds to a fold with acute angle. The degeneracy of each vertex under cyclic permutations 
of the links is indicated. 





Figure 10: Examples of 3c? octahedral foldings of an elementary hexagon and the corresponding vertices 
of Fig. | 
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Figure 11: Bending Energy vs rigidity for 3d folding from a transfer matrix calculation on strips of 
width L = 2,3,4,5 and 6. 



line corresponds to a fold with acute angle (cos -1 (—1/3) ~ 109°28', flip of z only). The degeneracy 
of each vertex under cyclic permutations of the links is also indicated. In Fig. ^ we display a few 
examples of vertices and the corresponding foldings in 3c? space. The 3c£ octahedral folding problem is 
thus precisely defined as a 96-vertex model with the vertices of Fig. ^. Although the entropy of the 3d 
problem cannot be evaluated exactly it can be determined by a numerical transfer matrix calculation to 
be q%d ~ 1.43(l).p5[ One can also derive various exact bounds on the entropy for 3d folding by relating 
it to dressed 3-colouring and 2d folding in an external staggered magnetic field. J|5| The best bounds 
derived in this way, expressed numerically to three decimal places, are 1.436 < q^d < 1.589. 

We have also examined the phase diagram of 3c? folding with bending energy. As for planar 
folding one finds a first-order phase transition to an ordered state at a critical value of the bending 
rigidity. The bending energy versus rigidity is shown in Fig. 11 for infinite strips of width L = 2, 3,4, 5 



and 6. One sees quite clearly the emergence of a non-zero latent heat as the system size increases. This 
is confirmed by a finite-size scaling analysis of the peak of the specific heat plot: it grows linearly with 
system size, as characteristic of a first-order transition. 

The planar and 3c? folding models described here have also been studied in the hexagon approx- 
imation of the cluster variation method. p6]] For planar folding this allows the incorporation of defects 
in the lattice and the study of the crossover from the pure Ising model to the pure folding problem. 
For 3c? folding these authors also add a symmetry-breaking field to the model and find a first-order 
transition to a flat phase for any value of the symmetry breaking field. 

In conclusion folding is a rich problem providing considerable insight into the phase structure 
of fluctuating membranes. 
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